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Confirmation of the nature of the absorber in IRAS 09104+4109 
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ABSTRACT 

We present the first long Suzaku observation of the hyperluminous infrared galaxy IRAS 
09104+4109 which is dominated by a Type 2 AGN. The infrared to X-ray SED indicates that 
the source is an obscured quasar with a Compton-thin absorber. However, the 3er hard X-ray 
detection of the source with the BeppoSAX PDS suggested a reflection-dominated, Compton- 
thick view. The high-energy detection was later found to be possibly contaminated by another 
Type 2 AGN, NGC 2785, which is only 17 arcmin away. Our new Suzaku observation offers 
simultaneous soft and hard X-ray coverage and excludes contamination from NGC 2785. We 
find that the hard X-ray component is not detected by the Suzaku HXD/PIN (effective en- 
ergy band 14-45 keV). Both reflection and transmission models have been tested on the latest 
Suzaku and Chandra data. The 0.5-10 keV spectrum can be well modelled by the two scenar- 
ios. In addition, our analysis implied that the absorption column required in both models is 
Nn ~ 5 x 10 23 cm -2 . Unless IRAS 09104+4109 is a "changing-look" quasar, we confirm 
that it is a Compton-thin AGN. Although the lack of detection of X-ray emission above 10 
keV seems to favour the transmission scenario, we found that the two models offer fairly sim- 
ilar flux predictions over the X-ray band below ~ 40 keV. We also found that the strong iron 
line shown in the Suzaku spectrum is in fact a blend of two emission lines, in which the 6.4 
keV one is mostly contributed from the AGN and the 6.7 keV from the hot cluster gas. This 
implies that the neutral line is perhaps caused by disc reflection, and the reflection-dominated 
model is more likely the explanation. The transmission model should not be completely ruled 
out, but a deeper hard X-ray spectrum observation is needed to discriminate between the two 
scenarios. 
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1 INTRODUCTION 

The hyperluminous infrared (IR) galaxy IRAS 09104+4109 (z = 
0.442 ) hosted in a cD galaxy in a rich cluster dKleinmann et al.l 
1988) is one of the most powerful (IR luminosity > 10 46 erg s _1 ) 
sources within z = 0.5. The active nucleus of IRAS 09104+4109 
was classified as an obscu red Type 2 quasar based on its opti- 
cal spectrum characteristic dKleinmann et al.lll988l : [Hines & Wills! 
1 19931 : iTran et al.l l 2000j). while the near IR spectrum gave a con- 
sistent classification Tsoiferetal.lll996h . The ASCA observation 
showed that the source is bright in X-rays, wi th a powerlaw spec - 
trum which is typica lly seen in most AGNs (Fabia n et all 1 1994) . 
Iffines & WilQ Il993l) analysed the VLA 1.4 and 5 GHz obser- 
vations and found a double-lobed radio source with straight jets 
extending along north-west and south-east of the nucleus. The 
ROSAT High Resolution Imager (HRI) observation revealed spa- 
tial extended X-ray emission with a central dip in the X-ray in- 
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tensity dFabian & Crawfordlll995l) . which has been confirmed to 
be one of two cavities in north-west a nd south-east of the core 
with a new 76 ks Chandra ob servation dHlavace k-Larrond o et al.l 
l2012h . lO'Sullivan et all J2012h presented new GMRT observations 
and cl aimed that these X-ray cavities coincide with the radio jets, 
while iHlavacek-Larrondo et al.l ([2012) indicated that the cavities 
are coincident with hotspots. The peaked X-ray surface brightness 
profil e resolved by ROSAT implied that the cluster has a s trong cool 
core dFabian & Crawfordlll995l;ICrawford & Vanderriestlll996l) . 

A strong emission line, probably associated with the reflection 
from helium-like iron, was discovered in the ASCA spectrum. The 
similar fluxes found in ASCA (extrapolated) and ROSAT observa- 
tions imply that the 0.1-10.0 keV X-ray emission is also likely con- 
tributed by the hot gas in the cluste r, though the strong iron line 
might come from a hidden nucleus. Franceschin i et al.l d2000l) re- 
ported a 3ct detection of a transmitted component above 10 keV by 
the BeppoSAX PDS, which indicated that the spectrum should be 
interpreted as a reflection continuum absorbed by a Compton-thick 
absorber. The equivalent width of the strong iron line given by the 
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9.1 ks Chandra observation dlwasawa et al.ll200lh is ~ 1-2 keV, 
which further supports the reflection-dominated scenario. However, 
the limited 0.1-10.0 keV X-ray spectra are not sufficient in dis- 
tinguishing the transm ission- and the reflection-dominated models. 
IPiconcellie t al. ( 2007) analysed the 14 ks XMM-Newton observa- 
tion and found no significant preference to these models. They in- 
dicated a possibility of "contamination" from nearby sources, as 
the field of view of the PDS instrument of BeppoSAX is large (1 
degree x 1 degree). Also, the 2-10 keV luminosity measured by 
the transmission-dominated model is consistent with the expected 
value on the basis of the bolometric luminosity. The absorber along 
the line of sight of the nucleus of IRAS 09104+4109 is either 
Compton-thin, or changed from Compton-thick to -thin between 
the BeppoSAX and XMM-Newton observations in five years. 

Vigna li et al. I J20T ]]) found that both Chandra and optical/mid- 
IR spectral analysis imply heavy, but not Compton-thick obscura- 
tion (i.e. that A^h > 10 24 cm -2 ). Moreover, they presented a 54- 
month Swift BAT map (15.0-30.0 keV band) which shows that the 
hard X-ray emission detected by BeppoSAX is probably associated 
with a z = 0.009 Type 2 AGN NGC 2785, which is only 17 arcmin 
away from IRAS 09104+4109. In order to investigate the ambigu- 
ous detection by BeppoSAX above 10 keV in depth, data from the 
hard X-ray band, where XMM-Newton and Chandra have no cov- 
erage, are necessary. 

The recent observation with Suzaku offers simultaneous X-ray 
data of 0.5-10.0 keV and 14.0-45.0 keV of the source. Given that 
the field of views of Suzaku X-ray Imaging Spectrometer (XIS) 
and PIN are 17.8 arcmin x 17.8 arcmin and 34 arcmin x 34 ar- 
cmin, respectively, emission from the nearby AGN NGC 2785 is 
likely included in the PIN detector if the pointing is aimed at IRAS 
09104+4109 (RA DEC = 09 13 45.49 40 56 28.2). The Suzaku PIN 
is a non-imaging instrument and emission from other sources can- 
not be simply excluded by region extraction. Therefore we shifted 
the pointing 6 arcmin west (new RA DEC = 09 13 21.49 40 56 
28.2) and assigned a 90° roll angle in our observation. The entire 
source still lies within the XIS detector after the change of point- 
ing, and contribution from NGC 2785 should be excluded from the 
field of view of PIN. 

Once the hard X-ray detection of IRAS 09104+4109 by the 
BeppoSAX PDS is also seen by the Suzaku PIN, it would confirm 
the reflection-dominated picture. Compton-thick AGN are impor- 
tant in understanding the evolution of suppermassive black holes. 
So far these ob jects are mostly found at low redshifts (z < 0.05; 
IComastrill2004l) . If IRAS 09104+4109 is Compton-thick, it will be 
one of the few distant Compton-thick AGN to have reliable X-ray 
spectroscopic data. Relatively-few such objects are known at the 
moment. Luminous Compton-thick AGN with well-studied multi- 
wavelength data are useful for understanding an important evolu- 
tionary phase in the growth of supermassive black holes. 

This paper presents analysis for the first Suzaku long observa- 
tion of IRAS 09104+4109 and offers conclusive results about the 
nature of the absorber in the source. It has been debated for years 
that the X-ray spectrum is transmission- or reflection-dominated, 
and we will discuss this in detail in section [4] Calculations in this 
paper were assumed a flat cosmology with Ho =71 km s ~~ 1 Mpc _ 1 , 
ttx = 0.73 and fi = 0.27. 
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Figure 1. The upper panel shows the PIN data (before background correc- 
tion) in black points and the total background in red points. The lower panel 
compares the Non-X-ray background and the earth-occulted background, 
which are shown in blue and green points, respectively. 



2 DATA REDUCTION 
2.1 Suzaku 

IRAS 09104+4109 was observed with Suzaku during 201 1 Novem- 
ber 18-20 resulting in a total of ~113 ks of good exposure time. 
The XIS was operated in the normal mode, and both 3x3 and 
5x5 editing modes were operated in all detectors (XIS0, XIS 1 and 
XIS3). As mentioned in the previous section, the pointing has been 
shifted by 6 arcmin. The source still lies perfectly within the XIS 
field of view, and no other point sources can be seen. Data have 
been reduced using the heasoft V6.11.1 software package with the 
latest calibration database following the standard procedure. The 
source spectrum was extracted using xselect from a circular re- 
gion with a radius of 110 arcsec. A larger circular region with a 
180-arcsec radius has been used for background estimation from 
a source-free region. Response files were produced by the script 
xisresp, which calls xisrmfgen and xissimarfgen to generate re- 
distribution matrix files (RMF) and ancillary response files (ARF), 
respectively. We combined spectra and response files of the two 
front-illuminated (FI) CCD XIS detectors (XIS0 and XIS3) using 
the script ADDASCASPEC in ftool. The spectral bins of all spec- 
tra were grouped for at least 20 counts. We use the FI spectrum 
over 0.5-10.0 keV, and the BI spectrum over 0.5-7.0 keV for fur- 
ther analysis, as the signal-to-noise ratio of the BI data is lower at 
high energies. 

The Hard X-ray Detector (HXD) was operated in XIS- 
nominal pointing mode. A non-X-ray background (NXB) and a 
cosmic X-ray background (CXB) should be combined to form the 
total HXD/PIN background spectrum. We obtained the NXB event 
file directly from the Suzaku Data Centre and extracted it using 
the background model D (the tuned model). The CXB was sim- 
ulated using the PIN response for flat emission distribution. The 
count rate of source before background correction is 0.340 ±0.002 
counts s -1 in the effective PIN energy band 14.0-45.0 keV (all the 
count rates of PIN will be given over this band if not specified), 
appearing below the total background level (0.343 ± 0.001 counts 
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s -1 ). As the count rates of the source and the background are at 
similar level (see also the upper panel in Fig. Q], we produced the 
earth-occulted background for a further check. The earth-occulted 
background is ideally consistent with the NXB, and it should re- 
place the NXB if discrepant. The lower panel in Fig. [TJ shows that 
the earth-occulted background (green data points) and the NXB 
(blue points) are in good agreement. The count rate of the earth- 
occulted background (0.344 ± 0.004 counts s _1 ) is slightly higher 
than that of the NXB (0.325 ± 0.001 counts s -1 ), confirming the 
non-detection in the PIN band. Thus we do not include the PIN 
spectrum in all our spectral fittings. Although the hard X-ray emis- 
sion from IRAS 09104+4109 is beyond detection in the HXD/PIN, 
we can estimate an upper limit of flux by including the systematic 
error of the NXB and CXB . We generated an ARF for the HXD/PIN 
using the hxdarfgen tool. The output ARF was coupled with the 
"hxdnominal" response file "ae_hxd_pinhxnomel 1_201 10601. rsp" 
to account for the off-axis correction. The resulting 14.0-45.0 keV 
fluxes of the NXB and CXB are around 2.16 x 1(T 10 ergs cm~ 2 
s _1 and 9.23 x 10~ 12 ergs cm~ 2 s _1 . Considering the systematic 
error of the NXB to be 3% (for observations with more than 10 ks 
exposure timeQ) and the lower limit of the total background flux 
can be calculated. The background-uncorrected 14.0-45.0 keV flux 
of the source is ~ 2.23 x 10 -10 ergs cm -2 s _1 , and we hence 
estimate the upper limit of the background-corrected source flux to 
be ~ 4.25 x 10 -12 ergs cm -2 s _1 . The observed BeppoSAX PDS 
11 ergs cm -2 s -1 , which is higher than 
12 ergs cm -2 s _1 in the same 



20-100 keV flux is - 
the PIN upper limit 
band. 



2.2 Chandra 
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The Chandra Advanced CCD Imaging Spectrometer (ACIS) ob- 
served IRAS 09104+4109 on 2009 January 6 for about 76 ks in 
very faint mode. We used tools in the latest CIAO V4.4 software 
package to reduce the data. The observation ends up with a good 
exposure of ~ 72 ks after background flare filtering. The nucleus 
of the cluster appears as a point source in the 0.5-7.0 keV (effec- 
tive energy band of the ACIS) image. Since IRAS 09104+4109 is 
a strong cool-core cluster showing a significant radial temperature 
gradient, background subtraction could be critical in some circum- 
stances. Thank to the excellent point spread function (PSF) of the 
Chandra ACIS, the nucleus spectrum can be extracted from a cir- 
cular region with a 1-arcsec radius. In this case, the background 
was extracted from the surrounding 2-4 arcsec annulus in order to 
avoid possible source emission outside the 1-arcsec-radius region 
and collect enough counts for the chi-statistics. The resulting spec- 
trum contains emission from the nucleus of the cluster, and has 
excluded the thermal contribution from nearby hot gas. We also 
generated a spectrum with a 55-arcsec radius source region and a 
55-110 arcsec annulus background region for comparison to our 
Suzaku data (regions used in the Suzaku observation are too large 
and would include edges of the chips of the ACIS). All the point 
sources resolved by Chandra within the background region have 
been removed for clean background extraction. Finally, we grouped 
both of our Chandra spectra with a minimum of 20 counts per bin 
as we did for the Suzaku data. 

The Chandra spectra obtained by different reduction strategies 
are shown in Fig. [2] There seems to be a drop between the < 2 
keV and the > 3 keV band in the nucleus spectrum, but it is not 
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Figure 2. The figure shows the spectra we extracted from both Suzaku and 
Chandra Observations. The data points have been mildly re-binned for clar- 
ity. The differences below ~ 1 keV between the Suzaku FI and BI XIS 
spectra are due to different calibrations. We ignore the Suzaku BI XIS data 
above 7 keV due to poor signal-to-noise ratio. 
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Figure 3. The figure shows the data/model ratio of the Chandra nucleus 
(magenta points) and 55-arcsec (blue points) spectra, which have been fitted 
with a simple powerlaw. It can be seen that the line energies are different. 



obvious in the spectrum extracted from the 55-arcsec radius region 
(hereafter the Chandra 55-arcsec spectrum). We also plotted both 
the Suzaku FI XIS and back-illuminated (BI) XIS spectra in Fig.[2] 
The shape of the Suzaku FI XIS spectrum is fairly similar to that 
of the Chandra 55-arcsec spectrum, which is expected because they 
were both extracted from a region including the nucleus and nearby 
hot gas. The flux difference between them is likely due to different 
background subtraction areas. A clear red-shifted iron line appears 
in all spectra. Interestingly, the iron line in the Chandra nucleus 
spectrum seems to peak at a different energy from that of the 55- 
arcsec and Suzaku spectra. This will be further investigated in the 
following section. 



3 DATA ANALYSIS AND COMPARISON 
3.1 Iron Line 

A st rong iron line with an equivalent width E W ~ 300-1200 eV 
(e.g. l Franceschini et al..2000l : llwasawa et alj200lHPiconcelli et al.l 
120071) has been reported in all previous X-ray studies of this source. 
The equivalent width is a model-dependent quantity, causing differ- 
ences in values quoted in previous literature. We first fitted a simple 
model composed of a powerlaw and a Gaussian line to the 3.0-7.0 



© 201 1 RAS, MNRAS 000.[TJ{8] 



4 C.-Y. Chiang et al. 



Table 1. The table lists fitting parameters that are bound in all data. JVh in this table shows the absorption column of ZPHABS. T is the photon index of the 
powerlaw component. kT\ and fcT 2 are the temperatures of the two thermal plasma components. Errors have been calculated at the 90 per cent confidence 
level. 



Model 


A^IO 22 cm" 2 ) 


r 


fcTi (keV) 


kT 2 (keV) 


Z/Zq 


X 2 /d.o.f. 


Reflection 
Transmission 


45 3+22.9 

~o c+23.0 
oa.o_ 21 g 


I .0+0.31 
1 ' o9 -0.19 
i OO+0.63 
1 - 28 -0.6S 


9 99 + I.I6 
z,z -0.43 
n 9C+I.I8 


a oc+1.64 
6 - 36 -0.90 

7 03+ 2 ' 97 
'• ua -1.67 


0.42t°:°06 
0.68±g;« 


1065.5/1152 
1072.4/1150 



keV band of our data. The rest line energy of the Suzaku spectra is 
6.631q:m keV, consistent with the value 6.64 ± 0.05 keV given by 
the Chandra 55-arcsec spectrum. Nevertheless, the rest line energy 
of the iron line in the Chandra nucleus spectrum is 6.42 ±0.06 keV 
instead of ~ 6.6 keV found in other spectra. We plot the data/model 
(fitted with a powerlaw only) ratio of the Chandra data in Fig. [3] 
in which the different peaks of the line can be seen. As the Chan- 
dra 55-arcsec and Suzaku spectra also include emission from the 
nucleus, the iron line shown in these spectra should be formed by 
two components, in which one of them is attributed to the nucleus 
and the other to the diffuse hot gas. In most literature to date, the 
iron line shown in the X-ray spectra of IRAS 09104+4109 has been 
treated as a single line instead of two. It is not surprising that the 
line energy and equivalent width of the iron line obtaine d from dif- 
ferent data sets are in disagreement. IVignali et al.l d201ll) also sug- 
gested the possibility that the iron emission line is originated by a 
blend of emission features from both the central AGN and cluster 
gas. 

We hereafter decomposed the iron line shown in the Chandra 
55-arcsec spectrum into two Gaussian components and re-fitted the 
data. The rest line energy of the second Gaussian component turned 
out to be 6.70^o'?o keV. The equivalent widths obtained for both 
iron line components are EW = 82I37 eV for the 6.42 keV line, 
and EW = 227±g 73 eV for the 6.70 keV line, respectively. As 
for Suzaku data, the equivalent widths are EW = 87±H eV for 
the 6.42 keV line and EW = 260±H 5 eV for the 6.70 keV line, 
which are consistent with the Chandra 55-arcsec spectrum. The 
equivalent width found in the Chandra nucleus spectrum is high 
(EW — 338^133 eV) possibly owing to a lower continuum level. 
Although in Fig.[3]it seems that the lines in both spectra are of simi- 
lar line widths, the line shown in the nucleus spectrum can be fitted 
by a narrow (line width < 0. 1 keV) Gaussian line and no further 
components are required. In the following analysis, we model the 
iron line in the 55-arcsec and Suzaku spectra with two decomposed 
line components. 

3.2 Spectral Fitting 

We started detailed spectral analysis by constructing a model to ex- 
plain the Chandra nucleus spectrum. The origin of the 6.42 keV Fe 
K line is generally attributed to reflected emission from cold mat- 
ter illuminated by high energy photons. Hence a neutral reflection 
continuum might be needed in the model. In addition, the struc- 
tures in the hard (> 3 keV) and soft (< 2 keV) energy band in the 
Chandra nucleus spectrum shown in Fig. [2] seems to be symmet- 
ric, implying that the < 2keV component can be the red-shifted 
Fe L complex caused by strong reflection. Ilwasawa et al.l (2001) 
also suggested this possibility in their analysis of the ~ 9 ks Chan- 
dra observation. We fit a model composed of an absorbed power- 
law component powe rlaw and a reflection continuum reflionx 
(Ross & Fabian 2005) in xspec to the Chandra nucleus spectrum. 
The photon index in the reflionx component has been set iden- 



Table 2. The table below lists the fitting parameters which have not been 
shown in Table[T] In this table Af pow , N me ^i, N me y^2, and JVrcf represent 
the normalisation of the powerlaw component, the two MEKAL components 
and the reflection component, respectively. F14-45 and -F20-8O are the 14- 
45 keV and 20-80 keV fluxes predicted by the model. All fluxes are shown 
in ergs cm -2 s~ 1 . * Since the predicted values for Suzaku XIS FI and BI 
spectra are fairly similar, we quoted that of the FI spectrum only. 



Parameter Suzaku Chandra Chandra 

FI & BI 55-arcsec nucleus 



Reflection 

Afpow 7.9t" 7 ° ° x 10~ 5 2A til x 

A^mcki S.ltti x 10- 4 1.3l°; 2 x 10- 4 

jV mek2 2.7l°; 4 x IO" 3 (fixed) 

N IC{ 3.3+^;° x IO" 5 3- 7 l?: 7 x 10- 5 

Fi4_ 45 *2.37 xl0~ 12 2.10 x IO" 12 3.20 X 10~ 12 

F20-80 *2.89 x IO" 12 2.55 x IO" 12 4.18 x 10~ 12 



Transmission 

TVpow 9-5t7 ? 4 5 x 10~ 5 l-6±l;f x 10- 4 

Nmokl 1 - 5 -l'.o x IO" 3 (1.4 ± 0.2) x 10~ 4 

Ar mok2 2 -°t{'° x 10~ 3 (fixed) 

Fi4_ 45 *1.89 x 10~ 12 1.68 x 10" 12 2.78 x 10~ 12 

F20-8O *3.12 x 10~ 12 2.76 x IO" 12 4.70 x 10~ 12 



tical with that in the powerlaw component, and allowed to vary 
in a range between 1.4 and 3.3 which is reasonable for the reflec- 
tion scenario. The ionisation parameter £ in reflionx was set to 
unity (£ = 1.0) to model the 6.4 keV Fe-K emission line from 
a neutral reflector with solar abun dance. The Galactic absorption 
column Nh = 1.8 x 10 20 cm" 2 jMurphv et ai1ll996h is also in- 
cluded and modelled by tbnewQ (Wilms, Juett, Schulz, Nowak, 
2012, in preparation) with the Wilms abundance. The absorption 
component in the source has been modelled using zphabs. Never- 
theless, the model cannot, even with a high iron abundance, explain 
the low-energy excess in the spectrum well. It seems more likely 
that the low-energy hump in the core-only spectrum is caused by 
photonionised gas. We added a thermal plasma component mekal 
(iMewe et alj|l985l ; iLiedahl et al.lll995l see also Kaastra 19923) 
which includes line emission from several elements and the Fe-L 
complex to the model and obtained a better fit. 

In order to extend the model to a version that can explain the 
Chandra 55-arcsec and Suzaku spectra, the 6.7 keV Fe XXV emis- 
sion line shown in these spectra should be modelled. If the line 
originates from reflection, it should come from part of a highly- 
ionised accretion disc which is usually close to the central black 
hole, where the illuminating source is nearby. By comparing the 
Chandra spectra extracted from different source regions, it is obvi- 

2 http://pulsar.sternwarte.uni-erlangen.de/wilms/research/tbabs/ 

3 http://www.sron.nl/files/HEA/SPEX/physics/megron.pdf 
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Figure 4. The upper set of figures shows the spectral fittings of the reflection and transmission models, while the lower set shows the decomposed components 
of the models. In the two upper plots, colours of each spectrum are expressed in the same way as in Fig. [2] In both lower figures, the black solid line stands 
the total model and coloured ones represent contributions from different components. The green and red dot-dash lines show the two MEKAL components of 
different temperatures. We plot the absorbed powerlaw as a blue solid line in both panels, and the reflection component has been shown as magenta dash line. 
The Galactic absorption has been taken out in these figures. 



ous that the 6.7 keV line is generated from the environment around 
the nucleus but not the central AGN. In addition, the X-ray emit- 
ting region of a typical AGN is expected to be small, and a spec- 
trum extracted from 2-4 arcsec, which corresponds to ~ 20-30 kpc, 
should contain no X-ray emission from the nucleus AGN. Hence 
the 6.7 keV is more likely scattered by the hot diffuse gas instead 
of reflected from the accretion disc. We then added another mekal 
component with a high gas temperature into the model to fit the 6.7 
keV line. The resulting full model can be expressed as: tbnew*( 

ZPHABS*POWERLAW + REFLIONX + MEKALi + MEKAL2). We set the 

normalisation of MEKAL2 to be for the Chandra nucleus spectrum, 
as this component is not required here. The Chandra and Suzaku 
observations differ by about two years, which is short for a quasar 
to evolve dramatically. We expected the Chandra and Suzaku data 
to be slightly different in flux only. The main variables in the model, 
that are, the absorption column Nn of zphabs, the photon index T, 
and the temperature kT and the metallicity Z/Zq in the mekal 



component, are bound in all data set. We assume that there is only 
an instrumental cross-calibration constant between the Chandra 
55-arcsec spectrum and Suzaku data, and bind normalisations of all 
model components together. As for the Chandra nucleus spectrum, 
normalisations of the powerlaw component, the reflionx compo- 
nent, and the MEKALi component are allowed to vary in a way dif- 
ferent from the other spectra. We summarised the fitting parameters 
of the reflection model in Table[T] and normalisations of parameters 
in Table [2] 

The reflection model we constructed is similar to the 
"absorbed powerlaw + mekal + reflection" model used in 
Franceschini et al. I l l2000l) . However, in their work an absorption 
column greater than 5 x 10 24 cm -2 is required, while in Table 
[T] it can be clearly seen that the absorption column required here 
is Ah ~ 4.5 x 10 23 cm" 2 , which is heavy but Compton-thin. 
The value is consistent with the number sugge sted by the trans- 
mission scenario (A H ~ 5 x 10 23 cm" 2 , Piconcelli eTai]|2007l : 
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Vignali et al]|201ll) . Since the HXD/PIN has no solid detection in 
the high-energy band, the transmission model is a plausible expla- 
nation as well. We constructed a transmission model by replacing 
the reflionx component with a Gaussian line using ZGAUSS in 
XSPEC. The constraint on the range of the photon index has been 
lifted. Again we bind the main parameters in all data as we did 
when fitting with a reflection model. We allow normalisations of 
the powerlaw component, the MEKALi component, and the ZGAUSS 
component in the Chandra nucleus spectrum to vary independently, 
while normalisations of these components in the other three spectra 
are bound together. The results have been again listed in Table[T]& 

HI 

The transmission scenario also suggests a Compton-thin ab- 
sorber with an absorption column of ~ 5.4 x 10 23 cm -2 , which is 
consistent with that obtained by the reflection model and previous 
studies. Most of the resulting fitting parameters of both models are 
close and consistent at the 90% confidence level. The photon index 
needed in each model is slightly different perhaps due to different 
broadband continuum (see also Fig.[4). The same reason may cause 
the differences in the temperatures of mekal components. Because 
the reflionx generates emission features in the soft X-ray band (~ 
1 keV) and the iron line band, temperatures required in the mekal 
components might be discrepant when the reflection continuum is 
replaced. The MEKALi component is probably contributed by the 
thermal emission in the nucleus (the 1-arcsec-radius circular re- 
gion), and the mekal 2 originated from the hot gas in the region 
outside the nucleus. The gas temperature of MEKALi we obtained 
by each model is below 4 keV, which is consistent with the results 
of lO' Sullivan et all d2012T) . Our fitting gives a gas temperature of 
~ 7 keV fo r mekal 2 , which also li es well in the 5-8 keV range 
indicated by lO' Sullivan et al.l d2012h . The values suggested by ei- 
ther the reflection or the transmission model are reasonable. As for 
the metallicity, values obtained by our models are close and only in 
mild disagreement. The best-fitted metalli city of the transmissio n 
model is closer to the value reported in O' Sull ivan et al. 1 d2012h . 
but that implied by the reflection model still lies within their 1 a 
uncertainties. 

The reflection and transmission models we tested in this data 
set result in comparable reduced x 2 ( see Table Q}. Even with a 
better signal-to-noise ratio, X-ray data cover across the 0.5-10.0 
keV band are not sufficient to distinguish these two interpretations. 
The lack of high-energy coverage makes model selection difficult, 
and the results are of no significant statistical difference. Fig. [4] 
shows the theoretical predictions from both models, and it can be 
clearly seen that they are fairly similar below 10 keV and start to 
divert above ~ 40 keV. The hard X-ray fluxes, especially the 20- 
80 keV band, predicted by these models are fairly close (see Table 
|2]l and both below the upper limit estimated in section [2] As the 
results of the two models statistically comparable, we will discuss 
in the following section the physical possibility of both models. 
The absorption-corrected rest frame 2-10 keV luminosities (cal- 
culated from the numbers shown in Table |3) of the nucleus are 
~ (1.3 - 2.3) x 10 44 erg s _1 and ~ (1.4 - 2.2) x 10 44 erg s"\ 
predicted by the reflection and transmission models, respectively. 
These values are slig htly higher than (1.2 — 1.3) x 10 44 erg s _1 
which was reported in lVignali et al.l ( l2011l) . but lower than the pre- 
dicted value of ~10 45 erg s " 1 obtained by t he sam e authors based 
upon IR observations from lLanzuisi et al.l d2009h . Using instead 
the m id-IR:X-ray correlation of local Seyferts from iGandhi et al.l 
(2009) results in an upper-limit to the X-ray power of ~ 5 x 10 
erg s -1 . This limit is consistent with the observations but sev- 
eral times higher, because of the fact that IR data on this source 



Table 3. The table lists the absorption-uncorrected rest from 2-10 keV 
flux *F2-\o, and the intrinsic 2-10 keV absorption-corrected flux F2— 10 ■ 
^thermal an d ^nucleus are the decomposed components of F2— 10 from 
thermal emission of the cluster gas and the nucleus. Fluxes are all shown in 
ergs cm -2 s — . 



Flux Suzaku XIS Chandra 55-arcsec Chandra nucleus 
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^2-10* 


2.00 x 10- 


■ 12 


1.76 x lO" 12 


1.04 x 10" 


■ 12 


-F2-IO 


1.74 x 10" 


■12 


1.54 x 10~ 12 


3.53 x 10" 


13 


-^thermal 


1.53 x 10" 


112 


1.36 x 10~ 12 


2.28 x 10" 


■14 


-^nucleus 


2.08 x 10" 


■ 13 


1.84 x 10~ 13 


3.31 x 10" 


13 








Transmission 






^2-10* 


2.14 x 10" 


■12 


1.89 x 10~ 12 


9.49 x 10" 


13 


F2-IO 


1.74 x 10" 


■12 


1.54 x 10~ 12 


3.53 x 10" 


■13 


-^thermal 


1.52 x 10" 


■12 


1.35 x 10~ 12 


3.77 x 10" 


14 


-^nucleus 


2.20 x 10" 


■13 


1.95 x 10~ 13 


3.15 x 10" 


■13 



are from the Spitzer mission, which cannot spatially resolve the 
nucleus emission from stellar ac tivity in typical ultra- and hyper- 
lumin ous infrared galaxies (see IVignali et alj|201 lL IGandhi et all 
l2009h . 



4 DISCUSSION 

4.1 Nature of the Absorber 

The main supporting evidence of the Compton-thick interpretation 
is the hard X-ray detection by the BeppoSAX PDS instrument. In 
our latest Suzaku observation, the pointing has been shifted to avoid 
possible contamination from the nearby Type II AGN NGC 2785. 
As a result the Suzaku HXD/PIN gave no solid detection but an 
estimated upper limit which is still lower than the measured Bep- 
poSAX flux, and the marginal 3a BeppoSAX detection is therefore 
questionable. The non-detection in th e hard X-ray band also im- 
plies that the reflection model used in Franceschi ni et alj fcOOOh is 
not necessarily required to explain the data. The hard X-ray flux 
obtained by the BeppoSAX observation was so high that the ex- 
trapolated transmission model failed to reach the flux, and the only 
possibility was to attribute the high-energy transmitted component 
to inverse Compton scattering. In order to explain the broadband 
X-ray data, the reflection continuum has to be heavily absorbed by 
extreme absorption column across the X-ray band below 10 keV 
and produce a Compto n hu mp to interpre t the ha rd X-ray emission. 
IPiconcelli etldl J2007l) and lVignali etail ( 1201 ll) obtained a similar 
high-energy flux prediction using a reflection continuum without 
absorber. They assumed that the absorber is of the same nature as 
the reflector, so using a reflection continuum automatically implies 
an optically-thick (r > 1) absorber. However, considering that the 
broad iron line displayed in the Suzaku XIS spectrum is a blend of 
two emission lines which can be caused by different mechanisms, 
the reflection scenario should not be ruled out. The narrow 6.4 keV 
Fe-K line could be reflected from the cold, optically-thick accre- 
tion disc in the central AGN. In this case, it is not controversial to 
model the data with a reflection continuum, as the optically-thick 
reflector could be part of the accretion disc and not necessarily act 
as an absorber. The reflection scenario is completely plausible here. 

In our analysis, the absorption columns required in the reflec- 
tion and transmission models are not extreme, and both scenarios 
give consistent values (~ 5 x 10 23 cm -2 ). The absorption col- 
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Figure 5. The figure shows the contour plot of column density of the ab- 
sorber TVjj against the photon index T. Contours are plotted at 67 (green), 
90 (red) and 99 (black) per cent levels. 

umn is correlated with the photon index, and we plot a contour 
of these parameters. It can been seen in Fig. [5] that the absorption 
column spans a range of Compton-thin values when the photon in- 
dex varies. In previous studies, the reflection model always links to 
a Compton-thick absorber, but we showed here that the reflection 
scenario does not necessarily imply the Compton-thick interpre- 
tation. As we mentioned before, the reflector needs not to be the 
absorber, and it is not surprising that we obtained a Compton-thin 
result using the reflionx model. If the BeppoSAX PDS detection 
is solid, the reflection model is the only interpretation capable of 
explaining the high-energy flux. Nevertheless, as the hard X-ray 
emission is weak, the reflection model can still explain the spec- 
trum well. The key to determine the nature of absorber is the flux 
above 10 keV, but not the model used to fit the data. It has been 
debated for a while that whether the absorber is Compton-thick or 
not. Now our analysis gives a conclusive result that the absorber 
in IRAS 09104+4109 is Compton-thin, provided that it is not a 
"changing-look" quasar. 

4.2 Reflection or Transmission above 10 keV? 

The two models we tested in this work gave statistically compa- 
rable results and similar fitting parameters. The predicted high- 
energy fluxes are both below the upper limit given by the 
HXD/PIN. We found that our predictio ns are also under th e up- 
per limit implied by the Swift BAT map dVignalietalJl201ll) . The 
3cr upper limit in the 15-30 keV band implied by the BAT map 
is ~ 1.9xl0~ 12 erg cm -2 s _1 , and the predicted fluxes in this 
bandare(1.2-1.8)xl0" 12 and (8.8-14.3)x 10" 13 erg cm -2 s" 1 for 
the reflection and transmission models, respectively. The reflection 
scenario has been ruled out in some studies due to its high flux 
prediction across the hard X-ray band. In our study, the reflection 
model works as well as the transmission model while no hard X- 
ray detection has been confirmed. The lack of high-energy emis- 
sion does not aut omatically imply that t he reflection model is out 
of consideration. IPiconcelli et al.l d2007l) indicated that the source 
does not display an iron line with a high enough EW expected for 
a truly reflection-dominated spectrum. Nonetheless, the reflection 
scenario a reflection spectrum does not necessarily come along with 
a strong iron line, as it is not the only feature produced by reflection. 



Iron Line 




Rest Energy (keV) 

Figure 6. The figure shows the data/model ratio of the Chandra 55-arcsec 
(blue open diamonds), the 2-55 arcsec (green triangles), and the 5-55 arcsec 
(black points) spectra, which have been fitted with a simple powerlaw. The 
red arrow points to the energy of the 6.4 keV line. 

For instance, the Compton hump and soft excess are signatures of 
reflection as well (though the soft X-ray band of IRAS 09104+4109 
is dominated by thermal emission from the cluster gas). In a case 
that the iron abundance is not high, the EW of the iron line can be 
low but the Compton hump is present. Thus the EW of the iron line 
should not be considered as an effective model selection tool. 

As mentioned in section 13.21 Fig. [4] clearly shows that both 
models give similar fluxes over the ~ 0.1-40.0 keV band. The re- 
flection model produces a Compton hump which peaks at ~ 40 
keV, while the transmission model predicts a steadily rising spec- 
tral shape above 20 keV. This implies that high-resolution X-ray 
spectroscopy in the hard-energy band may be required to tell which 
model provides a better interpretation to the data. If the spectrum 
is reflection-dominated, the shape of the Compton hump should be 
detected. As the HXD/PIN failed to detect emission above 10 keV, 
the only current space mission which is likely to achieve this is 
NuSTAR. 

Another way to examine the possibility of each model is to 
probe the origin of the neutral iron line. The 6.4 keV Fe Ka line 
could be fluorescent emission from (1) the accretion disc or (2) 
nearby low-ionised gas or molecular clouds. If the iron line is gen- 
erated from disc reflection, it must come from the nucleus of the 
source, where an AGN is accommodated. The line is likely re- 
flected from the edge of the cold outer disc as the torus. As for 
the latter poss ibility, a cooling flow can form clouds of atomic or 
molecular gas. lChurazov et ai1dl998l) indicated that the cold clouds 
illuminated by the X-r ay emission of th e hot gas would lead to a 6.4 
keV fluorescent line. lO' Sullivan et al.l d2012l) show optical emis- 
sion line filaments around the nucleus extended by a few arcsec, 
indicating that cold gas clouds are present. 

In order to examine the hypotheses mentioned above, we gen- 
erated Chandra spectra that exclude different sizes of central region 
and compared them with the Chandra 55-arcsec spectrum. We ex- 
tracted a spectrum which excludes the central 2-arcsec-radius circu- 
lar region, which means an extraction from a 2-55 arcsec annulus 
source region. Another spectrum excluding the central 5-arcsec- 
radius (5-55 arcsec annulus source region) was created using the 
same method. The background spectra were extracted using a 55- 
1 10 arcsec annulus, the same as the background region we used for 
the 55-arc spectrum extraction. We compare these spectra by plot- 
ting the line profiles in Fig. [6] and the red arrow indicates the 6.4 
keV line in the rest frame. It can be seen that there is a lack of strong 
emission at 6.4 keV rest frame in the 2-55 arcsec and the 5-55 arc- 
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sec spectra. We tested the transmission model on these new spectra 
and found that the 6.4 keV Gaussian line is not strongly required. 
The normalisation of the Gaussian component drops to a number 
close to zero, and the model works equally well if the Gaussian 
line is removed. This implies that most of the 6.4 keV emission is 
contributed by the very centre of the source, that is, the nucleus. 
The result cannot completely rule out the possibility that the 6.4 
keV line originates from the cold clouds, but shows that the neutral 
iron line is more likely caused by the AGN. We hereby consider 
the reflection scenario to be the better interpretation to explain the 
X-ray spectrum of this source. 



5 CONCLUSION 

We carried out a Suzaku observation which offers simultaneous 
soft and hard X-ray band monitoring. During the observation we 
avoided possible high-energy contamination from the nearby Type 
2 AGN NGC 2785. Our result does not confirm the 3<r detection by 
the BeppoSAX PDS instrument. By analysing both the latest Suzaku 
and Chandra observations with long exposures, we found that the 
broad iron line shown in these spectra is in fact a combination of 
two components. Both the reflection and transmission models have 
been tested on our data sets. Statistically the data are consistent 
with either of the interpretations and both models behave very sim- 
ilarly from the low-energy band to ~ 40 keV. Nevertheless, both 
models indicate a Compton-thin absorber in this source. Assum- 
ing that IRAS 09104+4109 does not change significantly within a 
short time scale, it is a Compton-thin AGN. Since the 6.4 keV Fe 
Kq seems to come from the central AGN, the reflection scenario is 
the better explanation to the source spectrum. However, the X-ray 
data below 10 keV give degenerate answers, and it is difficult to 
constrain the flux above 10 keV by current observations. We need 
high spatial resolution X-ray spectroscopy such as NuSTAR to de- 
tect the Compton hump predicted by the reflection model. 
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